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ABSTRACT 

Context. Emission from strongly-polar molecules could be a probe of physical conditions in diffuse molecular gas. 
Aims. We wish to provide basic information needed to interpret emission from molecules having higher dipole moments than CO, 
originating in diffuse clouds where the density is relatively low and the temperature and electron fraction are relatively high compared 
to dark clouds. 

Methods. Parameter studies in LVG models are used to show how the low-lying rotational transitions of common polar molecules 
HCO^, HCN and CS vary with number density, column density and electron fraction; with molecular properties such as the charge 
state and permanent dipole moment; and with observational details such as the transition that is observed. Physically-based models 
are used to check the parameter studies and provide a basis for relating the few extant observations. 

Results. Parameter studies of LVG radiative transfer models show that lines of polar molecules are uniformly brighter for ions, for 
lower J-values and for higher dipole moments. Excitation by electrons is more important for J= 1-0 lines and contributes rather less to 
the brightness of CS J=2-l lines. If abundances are like those seen in absorption, the HCO^ J=l-0 line will be the brightest line after 
CO, followed by HCN (1-0) and CS (2-1). Because of the very weak rotational excitation in diffuse clouds, emission brightnesses 
and molecular column densities retain a nearly-linear proportionality under fixed physical conditions, even when transitions are quite 
optically thick; this implies that changes in relative intensities among different species can be used to infer changes in their relative 
abundances. 
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1. Introduction 



The rotational excitation of CN in diffuse clouds was discussed 
as soo n as CN was identified in the diffuse ISM by i McKellaJ 
. (Il940l) . He calculated the "effective" or "rotational" temperature 

■ of interstellar space to be no more than 2.7 K and asked "if, 
indeed, the concept of such a temperature in a region with so 
low a density of both matter and radiation has any meaning." 
The efficacy of CN as an early probe of the cosmic microwave 

, background temperature was con firmed by its absence in several 
searches for mm- wave emission ('Thaddeusl 119721: ICrane et all 
.1989; Black & van Dishoeck, 1991), but gaining insight into the 

, structure of the absorbing gas requires more than upper limits on 

■ the excitation temperature of such a trace molecule. Moreover, 
[ actually detecting emission from various species among the 

many now known to exist in diffuse clouds may even permit 
imaging of the structure and chemistry of the host gas. 

CO emissi o n from diffuse clouds has long been known (see 
iKnapp & Jural d 19761) although several of the lines discussed 
there are telluric as noted in Liszt (2009)) and measurement of 
its excitation temperatures in absorption provides direct infor- 
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mation on the ambien t partial thermal pressure p/k ^ n(H2)TK 
= 1-10x10^ cm-^ K (ISmith et aI.L [1971 [Liszt & Lucasl[T998h . 
However, the mm-wave brightness of CO emission from diffuse 
gas is for the mo st part proportional only to the to t al CO column 
density N(CO) dGoIdreich & Kwanl [19741 iLiszli |2007|) . some- 
what leaving the determination of the physical properties of the 
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medium to other means, which for the most part remain to be 
explored. 

Emission fr om diffus e clou ds is accessible in 18cm 
OH emission dCrutcher, '1979^ and 9cm CH emission 
(Magnani et aI.L [l998. 2003), but the lines are weak, the spatial 
resolution is relatively poor and the e mission arises from states 
that are difficult to model or interpret ( Felenbok & Roue^ll996l : 
iLiszt & Lucasl[T996ll2002l) . Determination of the internal prop- 
erties of diffuse gas would be well-served by wider detection of 
mm-wave emission of molecules more strongly polar than CO, 
whose emission brightness might depend more directly on the 
number density and electron fraction in addition to the molecu- 
lar column density, and which could be observed with arcminute 
or better spatial resolution. Moreover, receivers have recently 
undergone drastic improvement in the 3mm band, so that the 
weakness of lines (see below) should not be such a barrier to 
future progress. 

Emission from strongly polar molecules has been seen 
in limited fashion in diffus e gas, in HCO^ (iLiszt & Lucas , 
Twl; i Lucas & Liszi [19961: [Uszl [19971 iKopp etalL [1996: 
Falgaron e et alll2006h . CS lPrdla et allll989h and HCN (iLiszt . 
11997 ). Where the HCO^ column density and optical depth are 
known in absorption, the weak HCO^ emission from opaque 
lines requires a small excitation temperature given by the line 
brightness, so that t he density ca n be disc ussed inde pendent of 
the column density dLucas & Lisz t. 1996; Liszt & Lu cas. 2000|). 
More generally, the density of colliding partners and the col- 
umn density of emitting molecules are not separable from ob- 
servations of a single transition, but surveys of the relative abun- 
dances of the various molecules in absorption might allow some 
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of this unfortunate degeneracy to be relieved. Observing the rel- 
ative intensities of various species could a means of deriving the 
ambient physical conditions if the relative abundances of those 
species are sufficiently well constrained by other means, and vice 
versa. 

The plan of this work is as follows. In Sect. 2 we lay out the 
basic elements of model calculations of rotational excitation by 
collisions with electrons and H2 in diffuse gas for species with 
different and varying permanent dipole moment and abundance. 
The first of these studies, described in Sect. 3, is a purely para- 
metric study using LVG radiative transfer (Goldreich & Kwan, 
1197 4) to limn the borders of the relevant parameter space. The 
second is a more physically-based study of excitation in models 
where the thermal and chemical balance and CO brightness are 
calculated in some detail and the results are relevant to interpre- 
tation of emission line observations of polar species when the 
only other or a priori information available on the host cloud is 
the strength of CO emission. The results of these studies are pre- 
sented along with a summary of the extant observations in Sect. 
4. The results of this work are summarized in Sect. 5. 



2. Models of excitation of polar molecules in diffuse 
clouds 

Molecule-bearing diffuse clouds are mar ked by elevated tem- 
peratures: <T| < > = 70-80 K for H2 (ISavage et al.L Il977t 
iRachford et al.L l2002h generally, while somewhat lower tem- 
peratures ar e found for regions with hi g her N(CO) - 10^- 
10'^ cm-2 (ISonnentruckeretal.L I2007L ISheffer et all l2007h . 
They are also marked by modest density n(H) = 30 - 500 cm~^ 
and high electron fraction n(e)/n(H) ^ 1 - 2 x 10"^, the latter 
corr esponding to the fr ee gas phase abundance of elemental car- 
bon (ISofia etal.Ll200l and a con tribution from cosmic -ray ion- 
ization of hvdrogen (^ LisztL 12003 ). Excitation calculations also 
must take into account that the J=l-0 rotational transition of a 
species having a permanent dipole moment /vDebye(lD = 10"'^ 
esu) will have t = 1 at column density per unit velocity dN/dV 
» 2 X lO'^/i"^ cm"^ (kms"')"', as is observed for many species. 

When the excitation is as weak as that which occurs in dif- 
fuse clouds for strongly-polar molecules, or said somewhat dif- 
ferently, so very sub-thermal (even in CO although it is much 
more strongly excited), emission line brightnesses are linearly 
proportional to column density even when the line optical depth 
is large. This occurs because, with so little collisional excitation, 
there is also little collisional de-excitation and most collisional 
excitations are followed by radiative de-excitation (Lyman-o' ra- 
diative transfer is an extreme example of the same phenomenon). 
Photons escape the medium even if they are scattered and re- 
absorbed many times (iGoldreich & Kwani 119741) . For CO the 
brightne ss is also insensitive to the ambient d ensity at fixed 
N(CO) (IGoldreich & KwanL 119741: iLisztl l2007h and the elec- 
tron fraction is unimportant because of the small dipole moment 
(0.1 ID). 

2.1. Excitation by electrons 

Excitation of strongly polar molecules by electrons has a 
strongly dipole character (AJ= +1) and the rates for such transi- 
tions are quite accurately represented in closed form for molecu- 
lar ions (Dickinson & Flower, 1981; Bhattacharvva et akl Il98lh 
and neutrals (Dickinson et al., 1977) having dipole moments 
fi <: 0.5D. The two formulations of the e-ion excitation rates 
are entirely equivalent although somewhat differently expressed 



(see also iNeufeld & Dalgarnd (Il989l) \ More accurate rates for 
e-HCO^ and e-HCN collisions including |AJ| <> 1 transi- 
tio ns have recently been calculated by Faure et al. (2007 aJ and 
Fau re et aP (l2007bl) . respectively, and for CS by Varambh ia et al.l 
(.20 101) and a comparison between calculations using them and 
the older rates in closed form is given below for HCO^ in Fig. 3. 



2.2. Excitation of strongly-polar molecules by neutral 
particles 

Excitation rates for Ha-HCO^ interactions were extensively 
recalculated by Flower (1999). Other species explictly dis- 
cussed here are C S (Green & Chapman, 1978) and HCN 
(iGreen & Thaddeusi [1974: .Monteiro & StutzkL ,1986). The de- 
pendence of these collisional rate constants on the dipole mo- 
ment of the target molecule is ignored here because they should 
be relatively weak compared to the other dependencies such as 
the electron-excitation rates and the spontaneous emission rates. 



2.3. Excitation of carbon monoxide 

Th e exci tation of C O in d i ffuse c louds was recently discussed 
by iLiszli (20063 and iLisztl (l2007b. Note that the very larg e H- 
CO excitation rates calculated bv lBalakrishnan et alJ (l2002h and 



extensively discussed by Liszt were subsequently repu- 

diated ([Shepler et al., 2007). Because of its low dipole moment 
CO is largely unaffected by the electron fraction and this is illus- 
trated in some of the parameter studies shown below where the 
molecular dipole moment is artificially varied to take advantage 
of the existence of analytic formulae for the electron excitation 
rates. For strongly sub-thermal excitation, the excitation temper- 
ature of the J = 1 -0 transition of CO is proportional to the ambient 
partial thermal pressure of H2, p/k = n(H2)TK, not to either Tk 
or n(H2) individually. 



2.4. Permanent dipole moments 

The permanent dipole moments used here are: CO, 0.1 12D; 
HCO+, 3.889 D; CS, 1.96 D; and HCN, 2.98 D. In our earlier 
work we used a value fi - 4.07 D for HCO^, which gave smaller 
HCO+ column densities by a factor (4.07/3.889)^ = 1.10 for a 
given integrated optical depth. 



2.5. Radiative transfer 

Strongly polar molecules typically become optically thick in 
their low-lying transitions for column densities around dN/dV ^ 
10'^ cm"^ (kms"')"', which is typical of actual observed abun- 
dances. Therefore it is necessary to treat the radiative transfer in 
some fashion. Sect. 3 presents a series of parameter stu dies in 
the usual LVG approximation (IGoldreich & Kwanlll974l) which 
amounts to solving the full, coupled, excitation rate equations 
in the escape probability approximation, taking q = (1 - e"^)/T 
as the photon escape probability for a transition whose optical 
depth is T and setting the line source function S,, = q Bv(Tcmb) 
-I- (1-q) Bv(Tex). In these expressions Tex is the transition exci- 
tation temperature and By is the Planck function. As noted by 
IGoldreich & KwanI (11974 . using the escape probability formu- 
lation transforms the rate equations so that they resemble the 
optically thin limit but with all the radiative rates (basically, the 
Einstein A-values) multiplied by q. 
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dN(HCO+)/dV [cm"2 / km s"^] dN(HCO+)/dV [cm-2 / km s"^] 

Fig. 1. Rotational excitation of HCO^ by H2 and electrons in diffuse gas at moderate density in the LVG approximation. The 
horizontal axis of all panels is the column density per unit velocity of HCO^. The top panels show the brightness temperature above 
the cosmic microwave background of the J=l-0 line , at left for fixed electron fraction Xg = ne/n(H) = 1.4 x 10""^ and varying 
hydrogen number density n(H) - 2n(H2), at right for varying electron fraction at n(H) - 256 cm"^. The lower panels show line 
brightness ratios in the same manner. At upper left the calculation for n(H) - 256 cm"^ is shown twice, for Tk = 15 (which lies 
below) and 35 K. Note that the lower panels in each column have a linear ordinate scale. 



Section 4 presents a more detailed calculation for small uni- 
form density gas spheres in which the radiative transfer is calcu- 
lated numerically in the microturbulent approximation. 



2.6. Observed relative abundances 

The abundances of HCO^, HCN and CS observed in absorption 
at mm- wavelengths are summarized in Fig. A.l of Appendix 
A. The calculations here were done with X(CS) = X(HCN) - 
X(HCO+) = 2 X 10-''. 



2. 7. Molecules not considered here 

The discussion here is limited to a few of the simplest species 
that are detected in diffuse gas. Failed searches for emission from 
CN ifi - 1 .45D) were noted in the Introduction but its intensity 
is severely diluted by hyperfine structure; HC N, whose abun- 
dance with respect to CN is fixed in diffuse gas (iLiszt & Lucasl 
l200lh . has a twice-larger dipole moment and is brighter in emis- 
sion, leaving little reason to c ontinue the old quest for CN. C2H 
is ubiquitous in diffuse gas (iLucas & Liszi 120001 : iGerin et all 



'201 ll), although its emission brightness will be diluted by hyper- 
fine structure. The dipole moment of C2H is 0.78 D. 

3. Parameter studies 

3.1. The HCO+ J=1-0 transition 

Figure 1 shows a composite of results for a calculation of HCO^ 
excitation, where the horizonal axis is the column density per 
unit velocity. As context, note that the HCO^ J=l-0 line has op- 
tical depth T = 1 for dN/dV - 1. 1 x lO'^ cm~^ (kmsy ') ' and 
many other mo l ecules like HCN jLiszt & Lucasl 200 Ih and CS 
(iLucas & Liszi l2002h increase in abundance relative to HCO^ 
for N(HCO+) ^ I - 2 X 10'^ cm'^. The largest values we 
have seen in absorption for in dividual kinema tic components is 
N(HCO+) a: 4-7 X IO12 cm-2 (iLucas & LisztL ll996) and typical 
line profile FWHM are 1-1.5 kms~'. A typical fractional abun- 
dance for HCO+ i s X(HCO+) = n(HCO+)/N(H2) = 2-3x10'' 
(iLiszt et al.L[2010l) . 

At the top in Fig. 1 the brightness of the J=l-0 line is plot- 
ted for a wide range of hydrogen number density at fixed Xe = 
njnu = 1.4x10""* and at right for varying Xe at n(H) = 256 cm"-'. 
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Note that n(H) = 2 n(H2) in the LVG calculations of Fig. 1-3 and 
the use of n(H) rather than n(H2) is for the sake of consistency 
with the more detailed calculations (see Sect. 4) in which n(H) 
is stipulated and the H-H2 equilibrium is actually calculated. 
Clearly the line brightnesses are linearly proportional to the col- 
umn density for both optically thick a nd thin emis sion as for 
sub-thermally excited CO (see Fig. 6 of lLisztl (12007 )). but they 
are also sensitive to the density of H2 and electrons. At left, the 
calculation for n(H) - 256 cm"^ has actually been repeated for 
Tk = 35 and 15 K to illustrate the weak temperature dependence 
of electron excitation of the J=l-0 transition. This simple calcu- 
lation accounts for the weak obs erved HCO ^ lines seen near ab- 
sorption line background targets dLucas & L iszt. 1996) with typ- 
ical values of the electron fraction at modest n(H) $ 300 cm"^. 
Note that t ^ 1 near the horizontal mid-point of each panel, so 
that T > 10 at the rightmost edge. Some deviation from linear- 
ity is seen at the highest densities and electron fractions where 
the excitation is stronger and more molecules are supported in 
higher levels of the rotation ladder. 

At top right in Fig. 1 it is seen that electron excitation ac- 
counts for some 75% of the line brightness under these nom- 
inal conditions. Electron fractions much higher than 10 are 
hard to maintain in denser gas given the free carbon abundance 
1.4 X 10 and the electron fraction in physically -based mod- 
els actually declines noticeably for n(H) > 100 cm"^, although 
carbon remains almost fully-ionized (see Sect. 4). 

In the lower panels we show the brightness of the higher ly- 
ing lines (the J=2-l line is often inaccessible near z = owing 
to atmospheric absorption). At moderate column density they 
are somewhat more dependent on the electron fraction than the 
number density (see at right) while the opposite is true at high 
density. At present, detecting emission from any but the J=l-0 
line of HCO^ seems unlikely for diffuse clouds. 

3.2. CSJ=2-1 

Figure 2 shows the brightnesses of the lowest three lines of CS, 
calculated using the close d form of the electr on-neutral molecule 
excitation rates (also see lDrdla et all (ll989h V As for HCO+ the 
CS lines decline strongly in brightness with each successive in- 
crease in J and although the HCO^ and CS J=l-0 lines will be 
about equally bright at the same column density, the predicted 
CS J=2-l lines are about a factor 6 weaker. About equal col- 
umn densities are typical for N(HCO^) ^ 10'^ cm"^ (see Fig. 
A.l). Weakness of the CS emission is due to a variety of fac- 
tors, CS being neutral and less strongly polar (1.96 vs. 3.89 D), 
but mostly the difference arises from comparison of J=2-l and 
J=l-0 transitions. Less of the line brightness is due to electrons 
in the CS J=2-l transition than for HCO^ J=l-0, perhaps half 
(top right), and very little at all for the CS J=3-2 line given the 
difficulty of moving population into higher J-levels with only 
AJ = +1 transitions. As shown in Fig. 2, the brightnesses of the 
higher-lying CS transitions are more sensitive to temperature be- 
cause more of the excitation is due to H2 and the rotation levels 
lie higher in energy so that purely thermodynamic considerations 
become more important. 

3.3. The importance of being polar 

To illustrate the influences at work in these excitation studies 
owing to the structure of the individual species, we performed 
another sort of parameter study which takes advantage of the 
availability of closed forms for the electron excitation rates. The 



excitation and line brightness are calculated for various transi- 
tions and molecules, but the permanent dipole moment is treated 
as a variable whereever it appears explicitly. In particular we 
consider the J=2-l line of CS and the 1-0 transitions of HCN 
and HCO+. 

Figure 3 shows the excitation temperatures (at left) and line 
brightness (at right) as the permanent dipole moment varies from 
0.1 to 8 D, for the default values employed in Fig. 1-2 (nn = 256 
cm~^, njnu = 1.4 x 10""*, Tk= 35 K), for hypothetical species 
whose energy level spacings etc. are the same as for HCO^, 
HCN and CS. Highlighted are the series of models at those 
values closest (actually very very close) to the true dipole mo- 
ments. For HCO^, results are also shown for the coiTect perma- 
ne nt dipole moment an d the more accurate recent excitation rates 
of iFaure et all (l2007al) . The HCO^ lines are somewhat brighter 
and more highly excited when the more detailed electron exci- 
tation scheme is employed; larger diff'erences would appear for 
the higher-lying lines. 

The behaviour seen in Fig. 3 is curious in some regards. 
Excitation temperatures decline monotonically with fi for opti- 
cally thin lines and then increase when both the dipole moment 
and optical depth grow large but none of this behavior is apparent 
in the line brightness: in all cases the lines are simply brighter for 
higher dipole moments. It is generally noted that electron excita- 
tion provides about the same degree of excitation for all species 
independent of fi because the excitation rates and spontaneous 
emission rates both vary as fr. However the current cases are 
not pure, more so for CS J=2-l, and electron excitation generally 
dominates only when the dipole moments are large even for the 
species which actually are more strongly polar. The decreasing 
excitation temperatures at small fi occur because the excitation is 
dominated by neutrals so that increasing /i mostly serves to cause 
higher but uncompensated spontaneous decay rates. Eventually, 
increasing the dipole moment raises the line optical depths to 
the point that radiative pumping props up the excitation temper- 
atures. 



4. Existing observations of I-ICO+, IHCN and CS 
emission 

To compare observations and datasets on something like a com- 
mon scale, we searched the literature for observations of polar 
species with accompanying '^CO measurements and these are 
shown in Fig. 4. In the following we denote the integrated bright- 
ness of the J=l-0 line by Wco, while those CS and HCN are 
W(CS) and W(HCN). 

4.1. HCO+ 

Shown in Fig. 4 are observations of H CO^ and '^CO J- 1-0 from 
the HCO^ a bsorption line survey of ILucas & Liszll (Il996h (for 
CO also see 'Liszt & Lucas (1998')) and emission line observa- 
tions of the cl ouds occultin g ^ Oph by Kopp et al. ( 1996) and 
iLisztl (Il997h . iFalgarone et al.1 (l2006l) only reported '^CO mea- 
surements to accompany their observations of HCO^ in Polaris. 
The general run of HCO^ brightnesses in the range W(HCO^) = 
0.01 - 0.1 K kms"' is very much in line with the model results 
in Fig 1 for n(H) = 256 cm"^ and N(HCO+) = 0.3 - 3 x lO'^ 
as typically observed. HCO^ appears in diffuse gas at a level 
W(HCO+) ^ Wco/100. 
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1012 10^3 10^2 10^3 

dN(CS)/dV [cm-2 / km s"^] dN(CS)/dV [cnT^ / km s"^] 

Fig. 2. Rotational excitation of CS by H2 and electrons in diffuse gas at moderate density, much as in Fig. 1 . From top down the 
panels show the brightness temperature of the CS J=l-0, }-2-l and J=3-2 transitions. At left the calculation for n(H) = 256 cm"^ is 
shown for 15, 35 and 55 K in all panels. The observed value toward BL Lac is indicated. 



4.2. HCN 

HCN emission was probably det ected at the position of peak 
HCO+ emission South of ( Oph dLiszti Il997h at about 1/3 the 
strength of HCO^; this datapoint appears in Fig. 4 at Wco - 6 
K km s"' . It is in line with expectations for approximately equal 
column densities of HCO^ and HCN as observed in absorption, 
see Fig. A.l. 

4.3. CS 

Observations of CS J=2-l were made bv iLisztl d 19971) at po- 
sitions around the sightline to f Oph. Where HCO^ emission 
was strongest 30' South of the star, with W(HCO+) = 0.09 K 
kmal^NaistixW^GO) = 6 K kms"', CS was undetected at a (3cr) 
level 5 times weaker. CS was also undetected toward the star, in 
line with expectations for the calculated excitation conditions in 
Figs. 1 and 2 with N(CS) « N(HCO+) as observed in absorption. 
More recently (to be published), we failed to detect CS emission 
at several other locations in the Southern portions of the ^ Oph 
cloud, including in directions showing 15 K CO lines. 

However, much stronger CS emission was reported by 

iDrdla et al. (1989) along several sightlines toward bright stars 
used as targets for optical absorption line studies. It is possible 
that the strong CS emission arises from material located behind 



the target stars, as is often the case for CO emission in such sur- 
veys (Liszt, 2009) but the accompanying HCO^ lines would be 
very bright indeed, and anomalously bright with respect to CO, 
if they were several times stronger than the observed CS. The 
N(CS)/N(HCO^) ratio may vary over a wider range than was 
apparent toward the set of background sources we surveyed in 
absorption, where N(CS)/N(HCO+) < 3 according to Fig. A.l. 

4.4. Physical models 

In practice, emission lines of CO may well be all that is available 
to infer cloud properties when HCO^ or other polar species are 
sought and that is the common meaure we applied to compare 
observations in Fig. 4. To bridge the gap between observation 
and the LVG parameter studies described in Sect. 3, we ex tended 
our earlier models of H2 and CO formation (lLiszi.l2007h to cal- 
culate the emergent line brightness of HCO^ and CS as well as 
CO. The models self-consi stently determine the thermal equi- 
librium ( Wolfire et al.i 119951) and equilibrium abundances of H2 
and CO as a function of radius in small spherical gas clots of 
given number and column density n(H) and N(H), under the as- 
sumptiorQ that CO forms from the ordinary thermal recombina- 
tion of a fixed relative abundance of HCO+, X(HCO+) = 2x10 ^ 



' This general methodology is known as "faith-based chemistry" 
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Fig. 3. Dependence of the rotational excitation upon permanent dipole moment yi (Debye) at n(H) = 2n(H2) = 256 cm , Tk= 
35 K and Xg = 1.4 x 10""^ as before. The horizontal axis is column density per unit velocity and vertical axes are line excitation 
temperatures (at left) and brightness temperatures above the cosmic background. The CS-like species at top has the same rotational 
energy spacing and H2-excitation rates as CS, but the dipole moment is varied wherever it appears analytically in expressions for 
the line optical depth, brightness etc. and electron-neutral molecule excitation rates. Panels in the middle row are for the strongest 
hyperfine component of the J=l-0 line of a neutral HCN analog, and the lowest two panels in each column pertain to the J=l-0 line 
of a molecular ion analog of HCO^ (compare with Fig. 1). In each row one series of calculations is shown by red dotted lines; by 
coincidence those use values very close to the actual dipole moments of CS (1.96 D), HCN (2.98 D) and HCO^ (3.89D) . The dotted 
blue cu rve in each of the lowest two panels uses the coiTect dipole moment and the recent electron excitation rates of iFaure et al.l 
(l2007ah 



The emergent line brightness is calculated for a m icroturbulen t 
medium of fixed linewidth. For instance. Fig. 6 of' Lisztl (l2007h 
shows how Wco N(CO) for 0.1 K kms ' < Wco 10 K 
kms^', 10"* cm-2 < N(CO) ^ 10'*" cm^^ observed) even 
though T ^ 1 at N(CO) ^ 10'^ cm^^. 

Shown in Fig. 4 are two series of models for HCO^ and for 
CS. Each curve is for a fixed central column density N(H) = 0.7 
or 2 xlO^' cm"^ and the internal number density of hydrogen 
increases to the right in marked steps of 2"^ beginning with n(H) 



= 128 cm at the left-most point. The relative abundances of 
CS and HCO+ were taken equal, 2 x 10 ^ The HCO+ and CS 
brightnesses are slightly lower at higher N(H) and fixed Wco 
because models with higher N(H) form CO more readily and 
have slightly lower electron fractions. Models with higher N(H) 
have brighter CS and HCO^ lines at the same number density 
(i.e. at their respective left-hand sides) but only because they 
have commensurately higher CS and HCO^ column densities. 
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Fig. 4. Integrated emission line brightnesses for HCO^ J=l-0 
(red boxes), CS J=2-l (green diamonds) and HCN J=l-0 (ma- 
genta parallelogram) observed in diffuse gas, plotted against 
the integrated brightness of CO J=l-0. Data sources are given 
in Sect. 4.1-4.3. Chained red and green curves are model re- 
sults as discussed in Sect. 4.4. For both HCO^ (red) and CS 
(green), calculations are shown for clouds with N(H) - 0.7 and 
2.0 xlO^' cm"^ with number density n(H) increasing from 
128 cm-3 at left in marked steps of 2^'^. X(HCO+) = X(CS) 
= 2 X 10"^ was assumed; the CS brightness increases linearly 
with increased abundance X(CS) at fixed Wco while that of 
HCO"^ declines very slightly, as discussed in Sect 4.4 of the text. 



The models reproduce the general run of observations 
of HCO+ but the calculated W(HCO+)/Wco ratios are ba- 
sically fixed by the HCO^ recombination chemistry and the 
empirically-determined HCO^ abundance and are very difficult 
to fine-tune. The W(HCO^)/Wco brightness ratio actually de- 
creases slightly when X(HCO^) increases because the CO abun- 
dance increases even faster. This is opposite to the linear, pro- 
portional variation of W(CS) with X(CS) that occurs because 
changing X(CS) does not affect Wco- 

The physically-based calculation confirms the insight from 
the LVG parameter studies that J=2-l lines of CS will be much 
dimmer than J=l-0 HCO+ when X(CS) = X(HCO+) but the 
two model curves for CS are almost coincident with each other 
in this parametrization (as is also nearly the case for HCO^). 
This means that the only way of effecting a strong change in 
the CS/HCO+ brightness ratio (at fixed Wco) will be to vary 
X(CS). Although CS and HCO^ were not measured along the 
same directions in Fig. 4 with strong CS emission, the existence 
of comparably bright HCO^ and CS lines at similar values of 
Wco strongly imphes that X(CS) » X(HCO+). 



5. Summary 

In Sect. 3 (see Figs. 1-3) we discussed a series of parameter 
studies based on LVG radiative transfer calculations for actual 
molecules (HCO^, HCN, CS) and a series of synthetic analogues 
with varying permanent dipole moments. In diffuse molecular 
gas at weak excitation, lines of polar molecules are uniformly 



brighter for ions, for lower J-values and for higher dipole mo- 
ments. Excitation by electrons is more important for ions and 
for J=l-0 lines and contributes rather less to the brightness of 
J=2-l lines. The fact that HCO^ is observed in its J=l-0 line and 
at about the same column density as HCN or CS should make it 
several times brighter in emission than HCN J=l-0 or CS J=2-l. 
Low levels of HCO+ J= 1-0 emission (Tb 0.01-0.03 K) should 
be ubiqitous in diffuse gas if (as observed) HCO^ is itself ubiq- 
uitous at column densities around N(HCO^) - 10'^ cm"^. Of 
course abundances play a role in whether a species is detectable 
and the relative abundances of HCO^ and CS or HCN differ 
(the CS/HCN ratio varies less, see Fig. A.l) and may vary by 
amounts which are also comparable to the differences in bright- 
ness induced by the other factors discussed here. 

Perhaps the most important point to take away from the sim- 
ulations is the linear proportionality between intensity and col- 
umn density that applies to all species and low-lying transitions 
when the excitation is as weak as it is in diffuse gas. This is 
true almost without regard for optical depth over the range of 
column densities that are typical of observed diffuse cloudfl . 
This proportionality means that model line brightnesses calcu- 
lated under given physical conditions can simply be scaled over 
a very wide range of molecular column density in many cases (a 
notable exception is when the brightnesses of HCO^ and CO are 
compared, see Sect. 4). But it also implies that when different 
species are observed along the same lines of sight (more pre- 
cisely, under the same physical conditions), their relative bright- 
nesses are set by their relative abundances. Variations in their 
relative brightnesses from position to position result from varia- 
tions in their relative abundances, because variations in physical 
conditions do not shift the pattern of relative brightnesses ap- 
preciably (compare Fig. 1 and 2). Observing, say, CS J=2-l and 
HCO^ J=l-0 lines of comparable brightness along the same line 
of sight really implies that N(CS)/N(HCO+) » 1, a situation 
that has never been observed in our absorption line work. 

We gathered and compared existing observations of mm- 
wave emission from polar species in Sect. 4. To put the differ- 
ent datasets on a command basis we plotted their brightneess 
against Wco, the integrated brightness of the J=l-0 CO line. 
HCO^ emission at a level 1-2% of Wco is typical. In line with 
the results of the parameter studies, CS emission is not seen in 
the gas around ( Oph, even when HCO^ emission is relatively 
strong and somewhat weaker HCN emission is probably de- 
tected at a level some three times weaker then HCO^. However, 
CS emission is detected in the directions of several early-type 
stars used for optical absorption line studies, with brightnesses 
that are quite comparable to those of HCO^ at similar Wco in 
other directions. This implies somewhat larger variations in the 
CS/HCO^ abundance ratio than were seen in absorption, and 
larger values of X(CS). 

To bridge the gap between these observations and the sim- 
ulations and parameter studies we calculated the brightnesses 
of HCO^ and CS lines in diffuse cloud models that previously 
were used to calculate the brightness of CO. In these models 
CO forms from the recombination of a fixed relative abundance 
X(HCO^) = 2 X 10 This chemical relationship couples the 
emergent brightnesses of the two species and largely fixes the 
ratio of their line brightnesses that otherwise is difficult to ad- 
just. The models reproduce the observed HCO^/CO brightness 
temperature ratio (1-2%). 



- This also applies to CO, although the discussion here is really con- 
cerned with more strongly polar species. 



8 



H. S. Liszt: The importance of being polar 



Plotting the observed brightness of polar molecules against 
Wco creates degeneracies in the parametric dependences of 
the brightness on other physical variables, such that curves of 
W(CS) vs. Wco nearly coincide even for cloud models that dif- 
fer appreciably in their hydrogen number and column densities 
(Fig. 4). That is, for a fixed set of relative abundances, the rela- 
tive brightnesses are fixed over a wide range of physical param- 
eter space. In that case, the relative brightnesses are determined 
chiefly by the various X(CS), X(HCN), etc, an insight we also 
gained from the preceding LVG-based parameter studies. 
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Appendix A: Relative abudances of l-iCO^, IHCN and 
CS 

Shown in Fig. A.l are the column densities of HCO^, HCN 
and CS as observed in absorption t oward extragalactic m m- 
wave continuum sources in the work of Lucas & Liszt (1996) for 
HCO+, Liszt & Lucas (2001.) for HCN and .Lucas & Liszt (20021) 
for CS. 



